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Выполнено исследование динамических свойств центробежных аппаратов для рас- 
пределения минеральных удобрений по поверхности поля применительно к обосно- 
ванию допусков на колебания параметров процесса. По допускаемой неравномер- 
ности распределения удобрений по полю рекомендуется вычислять допуски на рав- 
номерность подачи удобрений на распределяющий аппарат, математическое ожи- 
дание и среднее квадратическое отклонение угла бросания и дальностей метания. 
Вероятностная модель системы представлена последовательным соединением зве- 
на постоянного запаздывания и двух апериодических звеньев. Численным модели- 
рованием определены характеристики звеньев. 
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Введение. Функционирование сельскохозяйственных машин происходит в 
условиях непрерывно изменяющихся внешних условий. Влияние различных 
факторов, главным образом случайных в вероятностно-статистическом 
смысле, сказывается на показателях технологического процесса распреде- 
ления удобрений. Известно, что неравномерное внесение удобрений при- 
водит к снижению их эффективности. Снижение прибавки урожая и ухуд- 
шение его качества происходит за счет нелинейности функции урожая от 
дозы внесения удобрений и полегания растений на площадках с повышен- 
ными дозами внесения. Неодновременное созревание растений на площад- 
ках с различной дозой внесения удобрений приводит к ухудшению каче- 
ства урожая и увеличению потерь при уборке. 

Постановка задачи. Равномерность подачи удобрений на распределяю- 
щий аппарат зависит от свойств удобрений, частоты вращения дозирующих 
органов либо от скорости и шага планок подающего транспортера. Частота 
вращения распределяющего диска влияет на дальность метания частиц 
удобрений и закономерность их распределения вдоль радиусов зоны рассе- 
ва. Неравномерность вращения распределяющих дисков приводит к повы- 
шению неравномерности распределения удобрений по полю. 

Изменение фрикционных свойств удобрений приводит к нарушению 
симметричности распределения относительно линии движения машины и, 
как следствие, снижению равномерности рассева. 

Цель исследования — получить модели для исследования влияния 
условий функционирования на равномерность рассева удобрений по полю 
и для обоснования допусков на колебания параметров процесса. 
Методика исследований. Процесс распределения удобрений рассмотрен 
при следующих исходных допущениях и условиях: распределение удобре- 
ний на выходе из аппарата задано функцией плотности угла бросания Ка); 
средние скорости метания частиц и средние коэффициенты парусности не 
коррелированы с углом бросания а; дальность метания частиц р является 
случайной величиной с известной функцией плотности вероятностей Кр); 
диаметр диска по сравнению с размерами зоны рассева считается малой 
величиной. 
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Понятие «доза внесения удобрений» применяем к бесконечно ма- 
лым площадкам так, что доза д- определяется соотношением: 


т о 
где Дд-— масса удобрений, внесенных на площадку; ДЕ-— площадь выделен- 
ной площадки. 

Площадь, засеваемую аппаратом при отсутствии перемещения ма- 
шины, т.е. при И, = 0, называем зоной рас- 
сева [1,2]. При выполнении перечисленных 
ранее условий линии равной интенсивности 
в зоне рассева имеют почковидную форму 
(рис.1). Аппарат находится в начале коор- 
динат и выбрасывает удобрения в радиаль- 
ных направлениях. Выделим в зоне рассева 
площадку аЕ ограниченную двумя радиуса- 
ми, образующими угол да, и двумя окруж- 
ностями с радиусами ри р + др. Угол а от- 
считываем от линии движения машины, 
угол В- от поперечной оси координат. 

Вероятность попадания удобрений 
на площадку аЁР определяем как вероят- 
ность совмещения двух событий: попадания 
значения угла бросания в интервал аа и по- 
падания дальности в интервал ар. Так как 
дальность метания некоррелирована с углом а, то элемент вероятности ар 
определяем по формуле умножения вероятностей: 

аР = Ка) аа Кр) др. (2) 

Интенсивностью засева площадки называем предел отношения рас- 
хода удобрений &О, приходящегося на площадку, к ее площади аЁ при 
аЕ-› 0, т.е. 

ао 


Т= Нт —=. (3) 
аг> 0 аР 
Расход ЧО определяем умножением общего расхода О на элемент 
вероятности @Р. Тогда, учитывая, что аР = р ар да, 


ОИ (4) 


При движении машины со скоростью И, зона рассева перемещается 
с такой же скоростью относительно неподвижной площадки аЕ В этом слу- 
чае интенсивность внесения удобрений на площадке будет переменной. На 
рис.1 изображены линии уровней интенсивности / = соп$Ё, В центре зоны 
рассева интенсивность максимальна, а на ее краях — минимальна. 

Доза внесения д- на неподвижной площадке 4 после прохода раз- 
брасывателя определяется интегрированием: 








Рис.1. Схема определения дозы внесения удобрений 
горизонтальным аппаратом 





Г= 


п 
+= | АЕ, (5) 
[0] 


где & — время внесения удобрений на площадке. 

Для исследования динамических свойств аппарата необходимо вы- 
числить дозы внесения удобрений на площадках, расположенных вдоль 
оси У на равных расстояниях. Интегрирование необходимо выполнять в ре- 
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жиме единого времени. Единичный скачок фактора задается в едином вре- 
мени так, чтобы некоторые из площадок полностью засевались до скачка, 
некоторые в переходном режиме и остальные в новом установившемся ре- 
жиме. В едином времени задаются пределы интегрирования. Исследование 


выполнено в системе компьютерной математики (рис.2). 


м; 





Обозначения величин: 

Оа - подача удобрений на аппарат, кг/с; 

\Мт — скорость машины, м/с; 

Мо - математическое ожидание угла бросания, рад; 

со - среднее квадратическое отклонение угла бросания, рад; 
Мр - математическое ожидание дальности метания, м; 

ор - среднее квадратическое отклонение дальности метания, м; 


О (Ма, ва, Мр, ор, Х, УХ) — функция дозы внесения удобрений от координат площад- 


ки Х, У и параметров зоны рассева. 
Моделирование скачка расхода Оа 
1=0...200 \:=1-0,25+10 \5:=30 Оз:=1 Д:=1 
91(\):= Ос # У<\5 
((Оо+д)) # У>У\5 \1:=О1(\,) 


Рис.2. Моделирование скачка подачи удобрений 


Графическое изображение реакции системы на единичное ступен- 
чатое воздействие называется переходной характеристикой. Аналитиче- 
ское выражение А(Ё) переходной характеристики называется переходной 
функцией. Переходную характеристику (рис.3) получали решением уравне- 
ний (1)-(5), описывающих работу аппарата, и аппроксимированием пере- 
ходных характеристик. Для аппроксимирования переходных характеристик 


5-образного вида использовали выражение: 











- | ИВ о 
В (Е, = 1- ехй — | }} —- 
(ет ит т 
|: 
где Ти т - коэффициенты аппроксимации. 
1=0.. 100 У, = 1.0.5 4 10 


0.15 


9(0,08,10,2.5,1,\); 
9(0,08,10,2.5,10,\); 
9(0,03,10,2.5,15,%): 105 


0.1 
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Рис.3. Графики переходных процессов дозы 
при скачке подачи (Х= 1; 10и 15) 


226 


(6) 


Технические науки 








Продолжительность переходного процесса определяется размерами 
зоны рассева и скоростью движения машины. 

Переходные процессы изменения дозы внесения удобрений на ли- 
НИЯХ (рис.4), параллельных направлению движения агрегата, пересчитыва- 
ли в переходные характеристики путем вычитания от текущих значений 
дозы ее значения до начала переходного процесса и деления разности ор- 
динат дозы на величину скачка входного фактора. 


Моделирование скачка Мо: 
1 
Мос = 0.8 1=0.. 100  \5=50 \Мт = 2 {= — 
Ут 


у 
Ма) =| | Мао И < 


у 
(Моо+ АМ) ® > 
Ут 


1=0.. 140 \} = 1.0.5 +1 





0.08 


9м(0,0.3,10,2.5,1,\); 0.06 


9м(0,0.3,10,2.5,4,\); 


ве 004 
9м(0,0.3,10,2.5,8,\); 
9М(0,03,10,2.5,- 4,%); 2 
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Рис.4. Графики переходных процессов дозы 
при скачке математического ожидания угла бросания (Х=1; 4; 8и-4м) 


Результаты моделирования. Оценка точности аппроксимации показала, 
что при /7=1 начало переходного процесса аппроксимируется плохо. При 
т=2 вся кривая аппроксимируется удовлетворительно. 

Модель системы при 1=2 представляется последовательным со- 
единением звена постоянного запаздывания и двух апериодических 
звеньев. 

Значения коэффициентов передачи К, постоянной времени Ти вре- 
мени запаздывания т, приведены в табл.1-3. 

При скачке подачи положение и размеры зоны рассева не изменя- 
ются, поэтому переходные процессы одинаковы в сечениях, равноотстоя- 
щих от линии движения. 


Таблица 1 
Параметры переходных и передаточных функций при скачке Оа 





Таблица 2 
Параметры переходных и передаточных функций 
при скачке математического ожидания угла бросания М», 
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-5 -0,0450 13,9 1,11 


-0,0319 


Таблица 3 
Параметры переходных и передаточных функций при скачке ©, 





При скачке математического ожидания угла бросания зона рассева 
смещается в сторону от линии движения. Это приводит к увеличению дозы 
при положительных Х и к уменьшению её при отрицательных Х. Это видно 
в табл.2, где коэффициент передачи имеет разный знак. 

Скачок параметра М, приводит к изменению размеров зоны рассе- 
ва, но коэффициент передачи в этом случае примерно на порядок меньше, 
чем при скачке О, М, или б,. Скачок параметра б, вызывает слабые 


переходные процессы, поэтому не рассматривается в данной работе. 
Передаточная функция, полученная при этих допущениях, имеет вид: 
Г е"° Р 
= (7) 
(Тр+ 1)" 
Амплитудно-фазовую частотную характеристику системы (рис.5) полу- 
чали из передаточной функции заменой комплексной переменной рна /@ : 
раем 


’ (10 )= ею 1’ (8) 








1 15 2 ИЕ 
[^) 


Рис.5. Амплитудно-частотные характеристики аппарата 
(Мм=1; 2; 3; 0,5 м/с) 


Амплитудно-частотная характеристика определяет дисперсию дозы 
и дисперсию урожая в зависимости от дисперсий параметров 
Мб М, 365, идри 


[4 


ан в? — 
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где А=1+ 7? 0’; В=2 Тод. 

Последнее уравнение использовали для построения графиков функ- 
ции (0) и для обоснования допусков на параметры системы (см.рис.5). 

Если колебания одного из параметров имеют гармонический харак- 
тер с амплитудой А. и математическим ожиданием М. , то коэффициент 
вариации одной гармоники равен 

у, = Е: 1 
№12 м, 

Аналогично определим коэффициент вариации дозы: 


00. (10) 


А, 100 
ЖЕ р . 11 
м. (11) 
Разделив уравнение (11) на (10), получим: 
у А М, М 


ЕЕ) 
у А, М. м 





= (12) 
Ч 

Последнее соотношение использовано для обоснования допусков 
на неравномерность параметров процесса. Предположим, задана допусти- 
мая продольная неравномерность дозы у, = 10%. Требуется установить до- 


пуск на колебания подачи О. Кроме того, известно, что 


У 
Мо=М, Г, В, тогда Ту в т (13) 
т р 


2 


График зависимости (13) показан на рис.6. 





у@ 


Рис.6. Обоснование допуска на колебания подачи 


ФЕ Ё Т, = 1. 0.05+ 0.5. 003) 


Допуск на колебания подачи может превышать допуск на колеба- 
ния дозы внесения удобрений в десятки раз, если частота колебаний до- 
статочно высока (см.рис.6). Постоянная времени Т заметно влияет на ве- 
личину допуска только при повышении частоты колебаний. 

Выводы. При малой скорости движения машины колебания подачи с 
и 0,5с! полностью сглаживаются, а при движении со скоростью 3 м/с 
сглаживаются колебания с @ 3 с', полоса пропускания увеличивается. 
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Величина параметра Т передаточных функций сильно влияет на 
сглаживающую способность аппарата в области низких частот. На величи- 
ну параметра Т влияют продольные размеры сечений зоны рассева. В цен- 
тре полосы рассева этот размер меньше и параметр 7Т меньше. 

Скорость движения машины (см.рис.5) является важнейшим пара- 
метром амплитудно-частотной характеристики аппарата. С увеличением 
скорости увеличивается полоса пропускания частот. 

Увеличение размеров зоны рассева положительно влияет на 
способность аппарата фильтровать или сглаживать случайные колебания 
параметров процесса распределения удобрений по полю. 
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\. СНЕВМОУОГОМ, Т. ЧЁНАКНОТ\, Т. ЕУАЗНЕМКО 


ТНЕ ОУМАМТС РКОРЕВТТЕ$ ОЕ ТНЕ СЕМТВТЕУСАЕ ОЕ\МТСЕ 
РОК ТНЕ УРЕВЕТСТАЕ ЕМТЕКТМС ОЕ МТМЕКАЕ РЕКТТТ2ЕК$ 


Тре гезеагсН оЁ дупатс ргорег@ез оЁ {Не сепыТида! демсез Гог ЧзБиЧоп оЁ 
ттега! ГеИтег5 оп а Яе4 зиТасе адараБу о зиб$апНаНоп оГ Пе адпт- 
$1015 оп ЯисшаНоп ог {Пе ргосез$ рагатеег. ТАе адпт!5$!оп$ оп ргорогГ/опаЖу 
ог дампа ог {Ле Ге тег$ оп а 45 Би{та демсе, ап ауегаде дцадгаНс тс!та- 
Чоп ог ап апае ог {Игом/пд ап а гапде оЁ сазИпд аге гесоплтепаеч {о $е® {Пе 
поп-ипйогтКу ог 4 БиНоп ог Не Гей тег$ оп а Яеч. 

Тре ргобаЫе тоде! о! {Пе зу$ет 15 гергезетжеа Бу {Пе сопзеси ме соппесНоп 
оГа ШпК оЁ соп$КапЕ геагааНоп ап мо поп-рено!с пК$. ТПе йпК сВагаеп5- 
{с5 аге де егттеа Бу Пе питейс тодейпд. 
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